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Abstract 

A search for the standard model Higgs boson (H) decaying to bb when produced 
in association with weak vector bosons (V) is reported for the following modes: 
W(/^i/)H, W(ei/)H, Z{}ifi)H, Z(ee)H and Z(vv)H. The search is performed in a data 
sample corresponding to an integrated luminosity of 4.7 fb^^, recorded by the CMS 
detector in proton-proton collisions at the LHC with a center-of-mass energy of 7 TeV. 
No significant excess of events above the expectation from background is observed. 
Upper limits on the VH production cross section times the H — > bb branching ratio, 
with respect to the expectations for a standard model Higgs boson, are derived for a 
Higgs boson in the mass range 110-135 GeV. In this range, the observed 95% confi- 
dence level upper limits vary from 3.4 to 7.5 times the standard model prediction; the 
corresponding expected limits vary from 2.7 to 6.7. 
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1 Introduction 

The process by which the electroweak symmetry is broken in nature remains elusive. In the 
standard model (SM) HJEl the Higgs mechanism is considered as the explanation |l4]-0. The 
search for the Higgs boson is currently one of the most important endeavors of experimental 
particle physics. 

Direct searches by experiments at the Large Electron-Positron Collider (LEP) have set a 95% 
confidence level (CL) lower bound on the Higgs boson mass of wh > 114.4 GeV [.lOJ . Direct 
searches at the Tevatron exclude at 95% CL the 162-166 GeV mass range ITT) , and the ATLAS 
experiment at the Large Hadron Collider (LHC) excludes, also at 95% CL, the following three 
regions: mn ^ 145-206, 214-224, and 340-450 GeV l|T2tiT4ll . Measurements of the W boson 
and top quark masses at LEP and the Tevatron, combined with precision measurements of 
electroweak parameters at the Z pole, provide an indirect constraint of Wh < 158 GeV at 95% 
CL 115]. The most likely mass for the SM Higgs boson remains near the LEP limit, where the 
Higgs boson decays predominantly into bb. Experiments at the Tevatron have set 95% CL 
upper limits on the production cross section for a Higgs boson in this low-mass region. These 
limits range from approximately 4 to 10 times the standard model prediction, depending on 
the channels studied |IT6] - |22| . The observation of the H — > bb decay is essential in determining 
the nature of the Higgs boson. 

At the LHC the main SM Higgs boson production mechanism is gluon fusion, with a cross 
section of ~17pb for »jh = 120 GeV [23-39]. However, in this production mode, the detec- 
tion of the H — > bb decay is considered nearly impossible due to overwhelming dijet pro- 
duction expected from quantum-chromodynamic (QCD) interactions. The same holds true for 
the next most copious production mode, through vector-boson fusion, with a cross section of 
~1.3pb |,40-44J . Processes in which a low-mass Higgs boson is produced in association with a 
vector boson ||45l have cross sections of ~0.66 pb and ~0.36 pb for WH and ZH, respectively. 

In this Letter a search for the standard model Higgs boson in the pp — )■ VH production mode 
is presented, where V is either a W or a Z boson. The analysis is performed in the 110-135 GeV 
Higgs boson mass range, using a data sample corresponding to an integrated luminosity of 
4.7fb^^, collected in 2011 by the Compact Muon Solenoid (CMS) experiment at a center-of- 
mass energy of 7TeV. The following final states are included: W(f/i/)H, W(et/)H, Z(^/^)H, 
Z(ee)H and Z(vv)H, all with the Higgs boson decaying to bb. Backgrounds arise from pro- 
duction of W and Z bosons in association with jets (from all quark flavors), singly and pair- 
produced top quarks (tt), dibosons and QCD multijet processes. Simulated samples of signal 
and backgrounds are used to provide guidance in the optimization of the analysis as a func- 
tion of the Higgs boson mass. Control regions in data are selected to adjust the simulations 
and estimate the contribution of the main backgrounds in the signal region. Upper limits at 
the 95% CL on the pp — ?► VH production cross section are obtained for Higgs boson masses 
between 110-135 GeV. These limits are based on the observed event count and background es- 
timate in signal-enriched regions defined by the output discriminant of a boosted-decision-tree 
algorithm [46 1 (bdt analysis). As a cross-check, limits are also derived from the observed event 
count in the invariant mass distribution of H — > bb candidates (m(jj) analysis). 

2 CMS Detector and Simulations 

A detailed description of the CMS detector can be found elsewhere ||47| . The momenta of 
charged particles are measured using a silicon pixel and strip tracker that covers the pseudo- 
rapidity range |//| < 2.5 and is immersed in a 3.8 T solenoidal magnetic field. The pseudora- 
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pidity is defined as rj = — ln(tan(0/2)), where 9 is the polar angle of the trajectory of a particle 
with respect to the direction of the counterclockwise proton beam. Surrounding the tracker 
are a crystal electromagnetic calorimeter (ECAL) and a brass-scintillator hadron calorimeter 
(HCAL), both used to measure particle energy depositions and consisting of a barrel assem- 
bly and two endcaps. The ECAL and HCAL extend to a pseudorapidity range of |?/| < 3.0. 
A steel/ quartz-fiber Cherenkov forward detector (HE) extends the calorimetric coverage to 
\t]\ < 5.0. The outermost component of the CMS detector is the muon system consisting of 
gas detectors placed in the steel return yoke to measure the momentum of muons traversing 
through the detector. 

Simulated samples of signal and backgrounds are produced using various event generators, 
with the CMS detector response modeled with GEANT4 [48] . The Higgs boson signal samples 
are produced using POWHEG [49] interfaced with the HERWIG [50] event generator. The diboson 
samples are generated with pythia 6.4 [51]. The MadGraph 4.4 [52] generator is used for the 
W+jets, Z+jets, and tt samples. The single-top samples are produced with POWHEG and the 
QCD multijet samples with pythia. The default set of parton distribution functions (PDE) 
used to produce these samples is CTEQ6L1 |53|. The pythia parameters for the underlying 
event are set to the Z2 tune [54 1. 

During the period in which the data for this analysis was recorded, the LHC instantaneous 
luminosity reached up to 3.5 x 10^'^ cm^-^ s^^ and the average number of pp interactions per 
bunch crossing was approximately ten. Additional pp interactions overlapping with the event 
of interest in the same bunch crossing, denoted as pile-up events (PU), are therefore added in 
the simulated samples to represent the PU distribution measured in data. 

3 Triggers and Event Reconstruction 
3.1 Triggers 

Several triggers are used to collect events consistent with the signal hypothesis in all five chan- 
nels. Eor the WH channels the trigger paths consist of several single-lepton triggers with tight 
lepton identification. Leptons are also required to be isolated from other tracks and calorime- 
ter energy depositions to maintain an acceptable trigger rate. Eor the W(//v)H channel, the 
trigger thresholds for the muon transverse momentum, pj, are in the range of 17 to 40GeV. 
The higher thresholds are used for the periods of higher instantaneous luminosity. The com- 
bined trigger efficiency is ^90% for signal events that pass all offline requirements, described 
in Sectionjl] Eor the W(ev)H charmel, the electron pj threshold ranges from 17 to 30 GeV. The 
lower-threshold paths require two jets and a minimum requirement on an online estimate of 
the missing transverse energy, evaluated in the high level trigger algorithm as the modulus 
of the negative vector sum of the transverse momenta of all reconstructed jets identified by a 
particle-flow algorithm Il55l . These extra requirements help to maintain acceptable trigger rates 
during the periods of high instantaneous luminosity. The combined efficiency for these triggers 
for signal events that pass the final offline selection criteria is >95%. 

The Z(/^^)H channel uses the same single-muon triggers as the W(^i/)H channel. Eor the 
Z(ee)H channel, dielectron triggers with lower pj thresholds (17 and 8 GeV) and tight isolation 
requirements are used. These triggers are ?»99% efficient for all ZH signal events that pass 
the final offline selection criteria. Eor the Z(i/i/)H channel, a combination of four triggers is 
used. The first one requires missing transverse energy > 150 GeV and is used for the complete 
dataset. The other triggers use lower thresholds on the missing transverse energy (evaluated 
for these cases using all energy deposits in the calorimeter), but require the presence of jets. 
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One of these triggers requires missing transverse energy above 80 GeV and a central ( | | < 
2.4) jet with pj above 80 GeV, and the other two require the presence of two central jets with 
Pt > 20 GeV and missing transverse energy thresholds of 80 and 100 GeV, depending on the 
luminosity. The combined trigger efficiency for Z{vv)¥L signal events is ^98% with respect to 
the offline event reconstruction and selection, described below. 

3.2 Event reconstruction 

The reconstructed interaction vertex with the largest value of J^i pjj, where px, is the trans- 
verse momentum of the f-th track associated to the vertex, is selected as the primary event 
vertex. This vertex is used as the reference vertex for all relevant objects in the event, which 
are reconstructed with the particle-flow algorithm. The PU interactions affect jet momentum 
reconstruction, missing transverse energy reconstruction, lepton isolation and b-tagging effi- 
ciency. To mitigate these effects, a track-based algorithm that filters all charged hadrons that do 
not originate from the primary interaction is used. In addition, a calorimeter-based algorithm 
evaluates the energy density in the calorimeter from interactions not related to the primary 
vertex and subtracts it from reconstructed jets in the event [56 J . 

Jets are reconstructed from particle-flow objects [55J using the anti-A;^ clustering algorithm f57], 
as implemented in the fastjet package [58, 59] , using a distance parameter of 0.5. Each jet 
is required to be within |f/| < 2.5, to have at least two tracks associated to it, and to have 
electromagnetic and hadronic energy fractions of at least 1% of the total jet energy. Jet energy 
corrections, as a function of pseudorapidity and transverse energy of the jet, are applied [60[. 
The missing transverse energy vector is calculated offline as the negative of the vectorial sum 
of transverse momenta of all particle-flow objects identified in the event, and the magnitude of 
this vector is referred to as E^'®^ in the rest of this Letter. 

Electron reconstruction requires the matching of an energy cluster in the ECAL with a track in 
the silicon tracker [[6l1 . Identification criteria based on the ECAL shower shape, track-EC AL 
cluster matching, and consistency with the primary vertex are imposed. Additional require- 
ments are imposed to remove electrons produced by photon conversions. In this analysis, 
electrons are considered in the pseudorapidity range \ri\ < 2.5, excluding the 1.44 < \rj\ < 1.57 
transition region between the ECAL barrel and endcap. 

Muons are reconstructed using two algorithms [[62| : one in which tracks in the silicon tracker 
are matched to signals in the muon chambers, and another in which a global track fit is per- 
formed seeded by signals in the muon system. The muon candidates used in the analysis are 
required to be reconstructed successfully by both algorithms. Further identification criteria are 
imposed on the muon candidates to reduce the fraction of tracks misidentified as muons. These 
include the number of measurements in the tracker and the muon system, the fit quality of the 
muon track, and its consistency with the primary vertex. 

Charged leptons from W and Z boson decays are expected to be isolated from other activity 
in the event. For each lepton candidate, a cone is constructed around the track direction at the 
event vertex. The scalar sum of the transverse energy of each reconstructed particle compatible 
with the primary vertex and contained within the cone is calculated excluding the contribution 
from the lepton candidate itself. If this sum exceeds approximately 10% of the candidate pj the 
lepton is rejected; the exact requirement depends on the lepton rj, pj and flavor. 

The Combined Secondary Vertex (CSV) b-tagging algorithm [[63l is used to identify jets that are 
likely to arise from the hadronization of b quarks. This algorithm combines in an optimal way 
the information about track impact parameters and secondary vertices within jets in a likeli- 



4 



4 Event Selection 



hood discriminant to provide separation of b jets from jets originating from light quarks and 
gluons, and also from charm quarks. Several working points for the CSV output discriminant 
are used in the analysis, with different efficiencies and misidentification rates for b jets. For a 
CSV > 0.90 requirement the efficiencies to tag b quarks, c quarks, and light quarks, are approx- 
imately 50%, 6%, and 0.15%, respectively 1641 . The corresponding efficiencies for CSV > 0.50 
are approximately 72%, 23%, and 3%. 

All events from data and from the simulated samples are required to pass the same trigger 
and event reconstruction algorithms. Scale factors that account for the differences in the per- 
formance of these algorithms between data and simulations are computed and used in the 
analysis. 

4 Event Selection 

The background processes to VH production are vector-boson+jets, tt, single-top, dibosons 
{VV) and QCD multijet production. These overwhelm the signal by several orders of magni- 
tude. The event selection for the BDT analysis is based first on the kinematic reconstruction 
of the vector bosons and the Higgs boson decay into two b-tagged jets. Backgrounds are then 
substantially reduced by requiring a significant boost in the pj of the vector boson and the 
Higgs boson II65I , which tend to recoil away from each other with a large azimuthal opening 
angle, A^(V', H), between them. The boost requirements in the Z{£i)ii and WH analyses are 
Pj > 100 and pj > 150 GeV, respectively. The corresponding efficiencies for signal events are 
approximately 25% and 10%. For the Z(i/i/)H analysis the boost requirement is pj > 160 GeV. 

Candidate W — > £i/ decays are identified by requiring the presence of a single isolated lepton 
and additional missing transverse energy. Muons are required to have a pj above 20 GeV; 
the corresponding value for electrons is 30 GeV. For the W(ei/)H analysis, E^^®^ is required 
to be greater than 35 GeV to reduce contamination from QCD multijet processes. Candidate 
Z — )■ decays are reconstructed by combining isolated, oppositely charged pairs of electrons 
and muons and requiring the dilepton invariant mass to satisfy 75 GeV < tnu < 105 GeV. For 
Z candidates the electron pj is required to be greater than 20 GeV. The identification of Z — ?► i/v 
decays requires £^'®® > 160 GeV. The high threshold is dictated by the trigger and is consistent 
with a significant boost in the pj of the Z boson. The QCD multijet background is greatly 
reduced in this channel when requiring that the E^^^ does not originate from mismeasured 
jets. To that end, a A^(£^'®®,jet) > 0.5 radians requirement is applied on the azimuthal angle 
between the E^'*^ direction and the closest jet with pj > 20 GeV and \f]\ < 2.5. To reduce 
backgrounds from tt and WZ in the WH and Z(t/t/)H channels, events with additional isolated 
leptons, Nai, with pj > 20 GeV are rejected. 

The reconstruction of the H — > bb decay is made by requiring the presence of two central 
(I //I < 2.5) jets above a minimum pj threshold, and tagged by the CSV algorithm. If more 
than two such jets are found in the event, the pair of jets with the highest total dijet transverse 
momentum, piOj)/ is selected. After the b-tagging requirements are applied, the fraction of 
H — )■ bb candidates in signal events that contain the two b jets from the Higgs boson decay 
is near 100%. The background from V + jets and dibosons is reduced significantly through b 
tagging, and sub-processes where the two jets originate from genuine b quarks dominate the 
final selected data sample. 

The BDT analysis is implemented in the TMVA framework l'66'|. To better separate signal from 
background under different Higgs boson mass hypotheses, the BDT is trained separately at each 
mass value using simulated samples for signal and background that pass the event selection 
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described above. The final set of input variables is chosen by iterative optimization from a 
larger number of potentially discriminating variables. These include the dijet invariant mass 
m{ii), the dijet transverse momentum Pt(]])/ the separation in pseudorapidity between the two 
jets I A;/(jj) |, the transverse momentum of the vector boson pj{V), the maximum and minimum 
CSV values among the two jets, the azimuthal angle between the vector boson and the dijets 
A(p{V,li), and the number of additional central jets Naj. A signal region, where observed and 
expected events are counted, is identified in the BDT output distribution by optimizing a figure 
of merit that takes into account the level of systematic uncertainty on the expected background. 

Table [T] summarizes the selection criteria used in each of the five channels for both the BDT and 
the m(ji) analyses. For the cross-check ^(jj) analysis more stringent requirements are imposed 
on several of the variables used for the BDT selection. In addition, explicit requirements are 
made on A^(y, H) and on Naj. For each Higgs boson mass tested, nin, events are counted 
in a 30GeV window centered on the mean of the expected dijet mass peak. For the Z(££)H 
modes the dijet mass distribution is asymmetric and the window is centered 5 GeV lower than 
Wh, while for the WH and Z(i/i/)H modes the window is centered at nin. For these modes a 
higher pj boost requirement is made resulting in more coUimated b jets and a mass peak more 
symmetric around mn. For every channel, the ?«(jj) analysis was found to be about 10% less 
sensitive than the BDT analysis. 

Table 1: Event selection for the BDT analysis. Where applicable, the tighter requirements for 
the m(ji) analysis are listed in parenthesis. Entries marked "-" indicate that no requirement 
is made for that variable. The first two lines refer to the pj threshold on the leading (y'l) and 
sub-leading (;2) jets. CSVmax and CSVmin are the maximum and minimum b-tag requirements 
among the two jets. 



Variable 


W(£i/)H 


Z{U)H 


Z(t/t/)H 


pt(;i) 


> 30 GeV 


> 20 GeV 


> 80 GeV 


vAh) 


> 30 GeV 


> 20 GeV 


> 20 GeV 


PT(jj) 


> 150 (165) GeV 


> 100 GeV 


> 160 GeV 


Pt{V) 


> 150 (160) GeV 


> 100 GeV 




^miss 


> 35GeV[forW(ei/)H] 




> 160 GeV 


A</>(V,H) 


- (> 2.95) rad 


- (> 2.90) rad 


- (> 2.90) rad 


CSVmax 


> 0.40 (0.90) 


> 0.244 (0.90) 


> 0.50 (0.90) 


CSVmin 


> 0.40 


> 0.244 (0.50) 


> 0.50 


Nal 


= 




= 


Naj 


-(= 0) 


-«2) 




A(^(£^'s^jet) 






> 0.5 (1.5) rad 



5 Background Control Regions 

Appropriate control regions that are orthogonal to the signal region are identified in data and 
used to adjust the Monte Carlo simulation normalization for the most important background 
processes: W + jets and Z + jets (with light- and heavy-flavor jets), and tt. For each of the 
search channels and for each of these background processes, a control region is found such 
that its composition is enriched in that specific background process. The discrepancies be- 
tween the expected and observed yields in the data in these control regions are used to obtain 
a scale factor by which the normalizations of the simulations are adjusted. For each channel. 
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this procedure is performed simultaneously for all control regions. The background yields in 
the signal region from these sources are then estimated from the adjusted simulation samples. 
The uncertainties in the scale factor determination include a statistical uncertainty due to the 
finite size of the samples and an associated systematic uncertainty from the differences in the 
shapes of the distributions that could affect the estimate of the yields when extrapolating to 
the signal region. These systematic uncertainties are obtained by varying the control region 
selection criteria in order to select regions of phase space that are closer or further from the sig- 
nal region. The systematic uncertainty assigned covers the largest variation in the scale factor 
value found. The procedures applied in the construction of the control regions include revers- 
ing the b-tagging requirements to enhance W + jets and Z + jets with light-flavor jets, enforcing 
a tighter b-tagging requirement and requiring extra jets to enhance tt, and requiring low boost 
in order to enhance Vbb over tt. 

Consistent scale factors are found for each background process across the different channels. 
For tt, V + udscg, and Zbb production the scale factors are compatible with unity within their 
uncertainties (10-20%). For Wbb, the control region selected contains approximately 50% Wbb 
and single-top events, with the remainder being tt and W + udscg, which are well constrained 
by their own control regions. A choice is made to assign the observed excess of events in this 
region all to Wbb, leading to a scale factor of 2 for this background, while the estimate of single- 
top production is taken from the simulation. Reversing this assignment has a negligible effect 
on the final result of the analysis. The total uncertainty (excluding luminosity) assigned to 
the Wbb yield in the signal region is approximately 30%. This includes a 15% uncertainty on 
the extrapolation of the yield from the control region to the signal region, determined in data 
with the method outlined above. The systematic uncertainty assigned to the predicted yield 
for single-top production is 30%. The diboson background is taken from the simulation and a 
systematic uncertainty of 30% is assigned. 

For Z(i/t/)H the QCD multijet background in the signal region is estimated from data using 
control regions of high and low values of two uncorrelated variables with significant discrim- 
inating power towards such events. One is the angle between the missing energy vector and 
the closest jet in azimuth, A(^(£^'^*, jet), and the other is the sum of the CSV values of the two 
b-tagged jets. The signal region is at high values of both discriminants, while QCD multijet 
events populate regions with low values of either. The method predicts a very small contami- 
nation of 0.015 ± 0.008 for these background events, which is considered to be negligible. For 
all other search channels, after all selection criteria are applied, the QCD multijet backgrounds 
are also found to be negligible and not discussed in what follows. 

6 Yield Uncertainties 

Table |2] lists the uncertainties on the expected signal and background yields that enter in the 
limit calculation. 

The uncertainty in the CMS luminosity measurement for the dataset used in the analysis is esti- 
mated to be 4.5% [67]. Muon and electron trigger, reconstruction, and identification efficiencies 
are determined in data from samples of leptonic Z boson decays. The uncertainty on the yields 
due to the trigger efficiency is 2% per charged lepton and the uncertainty on the identification 
efficiency is also 2% per lepton. The parameters describing the Z(i/i/)H trigger efficiency turn- 
on curve have been varied within their statistical uncertainties and for different assumptions 
on the methodology to derive the efficiency. A yield uncertainty of 2% is estimated. 

The jet energy scale is varied within one standard deviation as a function of jet px and r]. The 
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efficiency of the analysis selection is recomputed to assess the variation in yield. Depending 
on the process, a 2-3% yield variation is found. The effect of the uncertainty on the jet energy 
resolution is evaluated by smearing the jet energies according to the measured uncertainty. 
Depending on the process, a 3-6% variation in yields due to this effect is obtained. An uncer- 
tainty of 3% is assigned to the yields of all processes in the WH and Z(vv)H modes due to the 
uncertainty related to the missing transverse energy estimate. 

Data-to-simulation b-tagging scale factors, measured in tt events, are applied consistently to 
jets in signal and background events. The measured uncertainties for the b-tagging scale factors 
are: 6% per b tag, 12% per charm tag and of 15% per mistagged jet (originating from gluons 
and light u, d, s quarks). These translate into yield uncertainties in the 3-15% range, depending 
on the channel and the specific process. 

The total VH signal cross section has been calculated to next-to-next-to-leading (NNLO) order 
accuracy, and the total theoretical uncertainty is 4% [39], including the effect of scale and PDF 
variations [68-72J. This analysis is performed in the boosted regime, and thus, potential dif- 
ferences in the pj spectrum of the V and H between data and Monte Carlo generators could 
introduce systematic effects in the signal acceptance and efficiency estimates. Two calculations 
are available that estimate the next-to-leading-order (NLO) electroweak |[73H75l and NNLO 
QCD 1 76 1 corrections to VH production in the boosted regime. The estimated effect from elec- 
troweak corrections for a boost of ~150GeV are 5% for ZH and 10% for WH. For the QCD 
correction, a 10% uncertainty is estimated for both ZH and WH, which includes effects due to 
additional jet activity from initial- and final-state radiation. The finite size of the signal Monte 
Carlo samples, after all selection criteria are applied, contributes 1-5% uncertainty across all 
channels. 

The uncertainty in the background yields that results from the estimates from data is in the 10- 
35% range. For the predictions obtained solely from simulation, as described in Section |5| an 
uncertainty of 30% (approximately the uncertainty on the measured cross section) is assigned 
for single-top. For the diboson backgrounds, a 30% cross section uncertainty is assumed. 

Table 2: Uncertainties in the signal and background yields due to the uncertainty in the sources 
listed. The ranges quoted are due to variations in mode, specific process, and Higgs boson 
mass hypothesis. See text for details. 



Source 


Range 


Luminosity 


4.5% 


Lepton efficiency and trigger (per lepton) 


3% 


Z(i/t/)H triggers 


2% 


Jet energy scale 


2-3% 


Jet energy resolution 


3-6% 


Missing transverse energy 


3% 


b-tagging 


3-15% 


Signal cross section (scale and PDF) 


4% 


Signal cross section (px boost, EWK/ QCD) 


5-10%/10% 


Signal Monte Carlo statistics 


1-5% 


Backgrounds (data estimate) 


10-35% 


Diboson and single-top (simulation estimate) 


30% 
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7 Results 

The primary physics result preser\ted in this Letter is an upper limit on the production of a 
standard model Higgs boson in association with a vector boson and decaying to a bb pair. 
Table |3] lists, for each Higgs boson mass hypothesis considered, the expected signal and back- 
ground yields in the signal region for the BDT analysis, together with the observed number of 
events. Table |3] also lists the requirements on the output of the BDT distributions that define 
the signal region. These distributions are shown in Fig. [T] for the nin = 115 GeV case, where 
data are overlaid with the predicted sample composition. The invariant dijet mass distribu- 
tion, combined for all channels, for events that pass the ?«(]]) analysis selection is shown in 
Fig. |2] The predicted number of background events are determined in data using the control 
regions described in Section |5| and from direct expectations from simulation for those back- 
grounds for which scale factors were not explicitly derived from control regions. Signal yields 
are determined from the simulations. The uncertainties include all sources listed in Section |6| 
except for luminosity. Total signal uncertainties are approximately 20%, and total background 
uncertainties are approximately in the 20 to 30% range. 

No significant excess of events is observed in any channel, and the results of all channels are 
combined to obtain 95% CL upper limits on the Higgs boson production cross section in the 
VH modes with H — > bb, relative to the standard model prediction. This is done separately 
for both the BDT and analyses for assumed Higgs boson masses in the 110-135 GeV range. 
The observed limits at each mass point, the median expected limits and the la and 2c bands 
are calculated using the modified frequentist method CLg l[77H79| . The inputs to the limit calcu- 
lation include the number of observed events (Npbs)/ and the signal and background estimates 
(Bexp)/ which are listed in Table |3] for the BDT analysis. The systematic and statistical uncer- 
tainties on the signal and background estimates, listed in Section |6| are treated as nuisance 
parameters in the limit calculations, with appropriate correlations taken into account. 

Table I summarizes, for the BDT and mij]) analyses, the expected and observed 95% CL upper 
limits on the product of the VH production cross section times the H — )■ bb branching ratio, 
with respect to the expectations for a standard model Higgs boson (a/ cgj^). The expected 
sensitivity of the BDT analysis is determined to be superior and it is considered to be the main 
result in this Letter. The BDT results are displayed in Fig.|3] 

8 Summary 

A search for the standard model Higgs boson decaying to bb when produced in association 
with weak vector bosons is reported for the following channels: W(^v)H, W(ei/)H, Z{mi)li, 
Z(ee)H and Z(t/i/)H. The search is performed in a data sample corresponding to an integrated 
luminosity of 4.7 fb^^. No significant excess of events above the expectation from background 
is observed. Upper limits on the VH production cross section times the H — )■ bb branching ra- 
tio, with respect to the expectations for a standard model Higgs boson, are derived for a Higgs 
boson in the mass range 110-135 GeV. In this range, the observed 95% confidence level upper 
limits vary from 3.4 to 7.5 times the standard model prediction; the corresponding expected 
limits vary from 2.7 to 6.7. This Letter reports the first upper limits from the LHC in these 
channels. 
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Figure 1: Distributions of the BDT output, for mH=115GeV, for each mode after all selection 
criteria are applied. The solid histograms for the backgrounds and the signal are summed cu- 
mulatively. The line histogram for signal is also shown superimposed. The data is represented 
by points with error bars. 
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Table 3: Predicted signal and background yields and observed niunber of events in data for 
the signal region defined by a BDT output value larger than the value listed. The uncertainty 

quoted is the total uncertainty, excluding luminosity. Results are given separately for each 
channel and Higgs boson mass hypothesis. W// and Z// denote W + udscg and Z + udscg, 
respectively. ST and YY denote single top and dibosons. 



W(ftv)H 



Process 


110 GeV 


115 GeV 


120 GeV 


125 GeV 


130 GeV 


135 GeV 


W// 


0.23 ± 0.14 


0.67 ± 0.29 


1.49 ±0.48 


0.39 ± 0.20 


1.48 ±0.48 


0.95 ± 0.38 


Wbb 


11.04 ±2.55 


7.78 ± 1.95 


8.32 ± 2.04 


4.50 ±1.30 


9.01 ±2.16 


6.89 ± 1.78 


Zbb 


0.84 ± 0.62 


0.84 ± 0.62 


1.29 ± 0.79 


0.84 ± 0.62 


1.29 ± 0.79 


1.29 ± 0.79 


tt 


1.66 ±0.59 


2.90 ±0.85 


2.90 ± 0.84 


1.31 ± 0.54 


2.79 ± 0.80 


1.71 ± 0.63 


ST 


1.32 ± 0.!52 


1.94 + 0.72 


2.58 ± 0.92 


1.74 + 0.66 


2.44 + 0.88 


1.59 + 0.60 


vv 


1.93 _ 0.66 


1.30 _ 0.46 


1.12 _ 0.40 


0.53 _ 0.21 


0.76 _ 0.29 


0.55 _ 0.21 


^exp 


17.02 ±2.83 


15.44 ±2.39 


17.70 ±2.60 


9.32 ±1.70 


17.78 ±2.65 


12.99 ±2.18 


Signal 


2.45 ±0.50 


2.06 ±0.42 


1.78 ±0.36 


1.08 ±0.22 


1.12 ±0.23 


0.75 ±0.15 


Nobs 


22 


23 


27 


15 


22 


13 


BDT 


0.21 


0.20 


0.20 


0.11 


0.21 


0.23 




W(ei/)H 


Process 


110 GeV 


115 GeV 


120 GeV 


125 GeV 


130 GeV 


135 GeV 


w;/ 


0.13 ±0.11 


0.51 ± 0.26 


0.37 ± 0.17 


0.23 ± 0.12 


0.28 ± 0.14 


0.10 ±0.11 


Wbb 


4.69 ± 1.06 


3.44 ±1.28 


3.72 ± 1.13 


3.53 ±1.10 


1.75 ±0.70 


2.08 ±0.76 


Zbb 


0.03 ±0.03 


0.03 ±0.03 










tt 


0.99 ± 0.46 


1.70 ±0.65 


2.31 ± 0.72 


2.07 ±0.71 


1.42 ±0.58 


1.17 ±0.51 


ST 


1.59 ± 0.59 


1.53 ±0.59 


1.75 ±0.67 


1.94 ±1.94 


1.51 ± 0.59 


1.33 ±0.52 


VV 


1.02 ±0.36 


0.63 ±0.24 


0.56 ±0.22 


0.45 ±0.18 


0.29 ±0.14 


0.25 ±0.12 


Sexp 


8.46 ±1.36 


7.84 ±1.59 


8.72 ±1.52 


8.21 ± 2.35 


5.25 ±1.10 


4.92 ±1.07 


Signal 


1.63 ±0.34 


1.39 ± 0.29 


1.20 ± 0.25 


1.04 ± 0.21 


0.76 ± 0.16 


0.61 ± 0.13 


Nobs 


9 


10 


10 


9 


8 


5 


BDT 


0.20 


0.13 


0.21 


0.09 


0.22 


0.24 





Process 


110 GeV 


115 GeV 


120 GeV 


125 GeV 


130 GeV 


135 GeV 


nf 


1.16 ±0.59 


0.95 ±0.52 


1.67 ±0.72 


0.62 ±0.42 


0.81 ± 0.48 


1.53 ±0.90 


Zbb 


4.85 ±1.48 


3.14 ±1.06 


7.05 ±1.98 


4.38 ±1.48 


5.67 ±1.79 


4.06 ±1.52 


tt 


0.64 ±0.22 


0.38 ±0.16 


1.05 ±0.32 


0.58 ±0.21 


1.19 ±0.32 


0.83 ±0.29 


VV 


0.92 ±0.35 


0.73 ±0.26 


1.01 ± 0.35 


0.55 ±0.20 


0.38 ±0.14 


0.15 ±0.06 


^exp 


7.57 ±1.64 


5.20 ±1.22 


10.78 ±2.16 


6.13 ±1.57 


8.05 ±1.88 


6.58 ±1.79 


Signal 


0.92 ± 0.17 


0.73 ± 0.13 


0.88 ± 0.16 


0.67 ±0.12 


0.59 ± 0.11 


0.43 ± 0.08 


Nobs 


7 


5 


6 


6 


11 


10 


BDT 


-0.207 


-0.195 


-0.246 


-0.221 


-0.313 


-0.243 



Z(ee)H 


Process 


110 GeV 


115 GeV 


120 GeV 


125 GeV 


130 GeV 


135 GeV 


Z// 


0.02 ± 0.02 


0.02 ±0.02 


0.20 ± 0.19 


0.32 ± 0.30 


0.36 ± 0.35 


0.02 ±0.02 


Zbb 


2.44 ±0.97 


2.51 ±0.98 


5.89 ±2.09 


5.48 ±2.04 


2.36 ±0.97 


3.44 ±1.19 


tt 


0.11 ± 0.08 


0.16 ±0.09 


0.38 ±0.17 


0.34 ±0.15 




0.12 ±0.09 


VV 


1.06 ±0.37 


1.07 ±0.38 


1.05 ±0.37 


0.92 ±0.33 


0.23 ±0.10 


0.46 ±0.19 


^exp 


3.63 ±1.05 


3.76 ±1.05 


7.52 ±2.14 


7.06 ±2.09 


2.95 ±1.03 


4.04 ±1.21 


Signal 


0.68 ± 0.13 


0.64 ± 0.12 


0.74 ± 0.14 


0.53 ± 0.10 


0.32 ± 0.06 


0.26 ±0.05 


Nobs 


2 


4 


4 


6 


5 


4 


BDT 


0.61 


0.63 


0.55 


0.59 


0.65 


0.67 




Z(vv)H 


Process 


110 GeV 


115 GeV 


120 GeV 


125 GeV 


130 GeV 


135 GeV 


w;/ 








0.89 ± 0.18 


1.53 ± 0.32 


1.53 ± 0.32 


Wbb 


4.46 ±0.99 


6.09 ±1.35 


6.12 ±1.35 


5.49 ±1.22 


3.23 ±0.71 


5.51 ± 1.22 


Z// 


1.27 ±0.24 


1.95 ±0.37 


1.20 ±0.23 


0.70 ±0.13 


0.66 ±0.13 


0.92 ±0.18 


Zbb 


5.74 ±1.42 


8.98 ±2.21 


6.47 ±1.30 


7.49 ±1.85 


8.77 ±1.76 


10.92 ±2.19 


tt 


1.04 ±0.12 


1.83 ±0.21 


1.96 ±0.23 


1.46 ±0.17 


1.19 ±0.14 


1.83 ±0.21 


ST 


0.61 ± 0.22 


0.85 ±0.31 


0.19 ±0.07 


0.27 ±0.10 


0.66 ±0.24 


0.53 ±0.19 


vv 


1.66 ±0.55 


1.64 ±0J4 


1.24 ±0.41 


0.56 ±0.18 


0.26 ±0.09 


0.41 ± 0.14 


Sexp 


14.78 ± 1.84 


21.34 ±2.70 


17.18 ±1.95 


16.86 ±2.24 


16.30 ±1.95 


21.65 ±2.55 


Signal 


1.82 ±0.33 


2.23 ± 0.40 


1.70 ±0.30 


1.64 ±0.29 


1.11 ± 0.20 


0.98 ±0.18 


Nobs 


15 


24 


20 


17 


16 


19 


BDT 


-0.18 


-0.17 


-0.15 


-0.20 


-0.22 


-0.25 
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Figure 2: Dijet invariant mass distribution, combined for all channels, for events that pass the 
m(jj) analysis selection. The solid histograms for the backgrounds and the signal are summed 
cumulatively. The line histogram for signal and for W backgrounds are also shown superim- 
posed. The data is represented by points with error bars. 



Table 4: Expected and observed 95% CL upper limits on the product of the VH production 
cross section times the H — > bb branching ratio, with respect to the expectations for a standard 
model Higgs boson. The primary results are those from the BDT analysis, the mij]) analysis is 
presented as a cross check. 



mH(GeV) 


110 


115 


120 


125 


130 


135 


BDT Exp. 
BDT Obs. 


2.74 
3.14 


3.12 
5.18 


3.56 
4.38 


4.27 
5.72 


5.28 
9.00 


6.74 
7.53 


7w(jj) Exp. 
OT(jj) Obs. 


3.05 
3.44 


3.21 
5.56 


4.14 
6.07 


4.74 
6.31 


6.42 
10.5 


7.67 
8.92 
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Figure 3: Expected and observed 95% CL upper limits on the product of the VH production 
cross section times the H — > bb branching ratio, with respect to the expectations for a standard 
model Higgs boson, for the bdt analysis. 
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